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We investigate a model for pattern formation in the presence of Galilean symmetry proposed by 
Matthews and Cox [Phys. Rev. E 62, R1473 (2000)], which has the form of coupled generalized 
Burgers and Ginzburg-Landau-type equations. With only the system size L as a parameter, we find 
distinct "small-L" and "large-L" regimes exhibiting clear differences in their dynamics and scaling 
behavior. The long-time statistically stationary state contains a single L-dependent front, stabilized 
globally by spatiotemporally chaotic dynamics localized away from the front. For sufficiently large 
domains, the transient dynamics include a state consisting of several viscous shock-like structures 
which coarsens gradually, before collapsing to a single front when one front absorbs the others. 

PACS numbers: 05.45.-a, 47.54.-r, 47.52. 02.30.Jr 



In the exploration of the rich and diverse range of spa- 
tiotemporal dynamics observed in nonlinear, nonequilib- 
riuni spatially extended systems, it has proved particu- 
larly fruitful to investigate comparatively simple model 
partial differential equations (PDEs) vifhose solutions 
capture the essential features of the phenomena under 
investigation. Thus the Burgers equation has been exten- 
sively studied for the evolution and statistics of shocks; 
the Ginzburg-Landau (GL) equation and its generaliza- 
tions describe the dynamics and stability of modulations 
of patterned states, while the Kuramoto-Sivashinsky and 
other models display spatiotemporal chaos (STC) [l[. In 
this paper we discuss a system describing the amplitude 
evolution for pattern formation with symmetry, which 
appears to combine features of several of these canonical 
systems and displays a surprising wealth of behaviors. 

We investigate the Matthews-Cox (MC) equations [2] 



At ^ A + AAxx - if A, 
It = fxx ~ \A\x 
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on a one-dimensional i-periodic domain, where A is com- 
plex, / is real, and fx = dxf = df/dX (similarly for 
the other derivatives). Equations ([I])-® were initially 
derived in the context of the Nikolaevskiy PDE 
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This equation, proposed originally to model seismic wave 
behavior in viscoelastic media Q, and subsequently ob- 
tained in other contexts [1, Q , appears to be a canoni- 
cal model for short-wave pattern formation with reflec- 
tion and Galilean symmetries. Unlike in more common 
pattern-forming contexts described at onset by the GL 
equation, the ©(e) stationary rolls in ^ are all unsta- 
ble for all £ > d, [H-01 ■ Instead, solutions of ^ ex- 
hibit spatiotemporal chaos with strong scale separation 



[si |9[, with coupling between the weakly unstable pat- 
tern at wave numbers fc « 1 and the neutrally stable 
long-wave mode with fc « 0. This suggests the Ansatz 
u{x, t) - e°'^A{X, r)e" + c.c. + e^/(X, T) + . . . for the 
envelopes A and / of the pattern and long-wave modes, 
respectively, where X = ex, T = e^t. Matthews and Cox 
showed that the asymptotically self-consistent scal- 
ing as e is ai = 3/2, /3 = 2, and hence derived 
©^(El) f^rom the Nikolaevskiy PDE as the leading-order 
modulation equations. While the scaling behavior on the 
attractor of ([3]) may be insufhciently described by this 
Ansatz ^9] , the MC equations deserve study in their own 
right as generic amplitude equations for pattern forma- 
tion with these symmetries [2]. Since ([2]) preserves the 
spatial mean of /, by Galilean invariance we may assume 
/ to have mean zero. 

Chaos- stabilized fronts: In describing properties of 
the MC equations (U)-©, we emphasize the dynamics 
of the large-scale mode /, since the pattern amplitude 
A appears to be driven by /. We note that several as- 
pects of the behavior for relatively small L have been 
previously described by Sakaguchi and Tanaka [lO^. 

The snapshot of a solution for domain size L = 51.27r 
shown in Fig. [T] is typical of the statistically stationary 
behavior for "small" domains. The overall structure of / 
resembles a perturbed viscous shock, with / decreasing 
essentially linearly within the "front" region. Simultane- 
ously, 1^1 vanishes in the center of the front; Sakaguchi 
and Tanaka hence call this an "amplitude death" state 
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FIG. 1: Snapshot at a fixed time Ti = 28000 of f{X,Ti) 
(thick blue line) and |^(X, Ti)| (thin black line) for a solution 
of Q-dll) with L = 51.27r 160.8. 
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FIG. 2: Space-time plots of long-time solutions (a) f{X,T) 
and (b) \A{X, T)\ of the MC equations for L = 51.27r. 

[Toj . The time evolution of a typical solution shown in 
Fig. [2] clearly shows the invariance of the front structure 
in / and the suppression of the roll amplitude A within 
the front region. On the remainder of the domain, chaotic 
fluctuations in / (on 0(1) time scales) are superimposed 
on the approximately linear average positive slope, and 
correlated (in space and time) with chaos in A. This 
coexistence of an ordered front (amplitude death state) 
and a spatially localized spatiotemporally chaotic region 
is robust on all domains L > Lq large enough to sustain 
the front ,10,] . 

The space-time plot indicates that the overall viscous 
shock profile in / is nonstationary, but maintains its 
shape up to small fluctuations; that is, short-time av- 
erages [12j {f{X, ■)}r are invariant up to translation. De- 
noting the averaged profile by g{X) — g''^\X), where we 
center the front so g{L/2) — 0, gx{L/2) < 0, and defin- 
ing the front displacement s{T) so that the instantaneous 
front position is Xs{T) = L/2 + s{T), we may decompose 
the large-scale mode as f{X, T) = g{X - s{T)) + f{X - 
s{T),T), where / denotes fiuctuations about the mean 
profile. The unsteady dynamics in A and / are then es- 
sentially confined to the region where gx > and to the 
vicinity of the local extrema of g (see Fig.[3l^c) below). 

To help clarify this unusual behavior, we observe that 
the equation ^ for the large-scale mode / has the form 
of a conservation law [13] , 

fT = -Jx, where J=-fx + \A\''. (4) 

Taking long-time averages, {J)x = —{It) = 0, so in sta- 
tistical equilibrium, the timc-avcragcd fiux J is uniform 
in X, (J) = {—fx + = 7- Integrating over the 

domain and using periodicity, we find 

l = l{L) = ^j\\A{X,-)\')dX> 0. (5) 

Now for a stationary amplitude death domain, where 
(|Ap) = 0, we have 7 = —{fx)', this confirms that am- 
plitude death can occur only where / is decreasing on 
average. (In fact, due to the apparent separation of time 
scales between the rapid fiuctuations / (on 0(1) times) 
and the slow overall drift of the mean profile g{X), short- 
time averaging seems sufficient to conclude {fT)T ~ 0, so 
that the mean mid- front slope is gx{L/2) = —7.) By 



^ the (averaged) slope —7 of the large-scale mode / in 
the center of the amplitude death state thus seems to be 
globally determined, being balanced by the mean-square 
amplitude of the pattern mode A due to chaotic dynam- 
ics concentrated outside the front region. 

Note that also 7 — {\A\'^) wherever {fx) = 0, relating 
the front slope to the fluctuations in A at extrema of the 
averaged proflle. Indeed, in the absence of A, f satisfles 
a heat equation, and thus by (HI flows away from lo- 
cal maxima and towards local minima, destabilizing the 
front; the added forcing term in ([2]) when ^ in- 

creases the flux J so as to maintain the averaged over- 
all flux balance, thereby stabilizing the (averaged) local 
extrema. The stabilization mechanism in the MC equa- 
tions thus appears to act globally [l^, with the chaotic 
dynamics being essential to sustaining the front (and am- 
plitude death state); we denote the observed structures 
"chaos-stabilized fronts". 

(We remark that ([1]) does not contain the usual sta- 
bilizing GL cubic term, permitting (UJ-© to support 
a family of exponentially growing solutions A{X, T) — 
Aqc^ , f{X,T) = 0; that is, the MC equations do not 
have a bounded global attractor. However, these grow- 
ing solutions are dynamically unstable, in the sense that 
they are overtaken by faster-growing spatially varying 
perturbations and in our numerical simulations we 
have not observed such solutions.) 

To investigate the behavior of the MC equations 
on spatially periodic domains systematically, we have 
numerically integrated ([I])-© using a pseudospectral 
method in space and an exponential time differencing 
(ETDRK4) scheme with step size H = 0.02. The do- 
main length L, the only free parameter in the system, 
was chosen to be L = 27r x 64m/10 for integers m ranging 
from 2 to 64; correspondingly, we used between 2^ and 
2^^ Fourier modes. In computing time averages, we inte- 
grated until the system reached a statistically stationary 
single-front state, and then averaged over 10^-10"* snap- 
shots separated typically by time intervals AT = 10. All 
averaging was done within the frame of reference of the 
front; that is, we first determined the front displacement 
s(T) and used it to align A and / so that the front was 
centered at X — L/2. In particular, the mean profiles 
were computed by g{X) = {f{X + s(-), •)). 

Averaged profiles: As seen in Fig. [31 the averages of / 
and \A\ are, respectively, odd and even about X — L/2, 
recovering the refiection symmetry of the underlying 
PDEs (HI)-©. More interestingly, though, the time- 
averaged profiles g{X) for L > Lo depend strongly on 
L, with the behavior falling into three distinct regimes: 

For relatively "small" domains, Lq < L < Li 220, 
the scaled profiles in Fig. [31[a) approximately coincide, 
indicating a scaling form for g: for some fundamental 
shape function G, periodic on [0,1], we have g{X) k, 
LG{Y) (with Y = X/L). In this "small-L" regime the 
scaling relation is highly accurate within the front (but 
is weakly violated outside it: the slope a = gx{0) m the 
active region increases slowly with L; see Fig. IHa)); in 
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FIG. 3: (a), (b) Two representations of the long-time- 
averaged (centered) profile g(X) = {f{X + s(-), ■)) of the 
large-scale mode / for various L: (a) scaled horizontally and 
vertically, L'^giLY) for Y = X/L e [0, 1]; (b) unsealed, half 
of the (odd) profile, g{X) for X e [0,L/2]. (c) Centered 
time-averaged pattern amplitude (|j4(I/Y, •)|'^). Domain sizes 
are ("small": red) L = 25.67r, 38.47r, 51.27r, 64tt; ("inter- 
mediate": green, dashed lines) L = 76.87r, 89.67r, 102. 47r, 
1287r, 153. Btt; and ("large": blue) L = 204.87r, 2567r, 307.27r, 
358.47r, 409. 67r, 460.87r, 5127r and 819. 27r. Averages were 
taken over time periods T = 1 x 10^ (small, intermediate) 
and 2 x 10^ (large), with AT — 10 between snapshots. 



particular, the midpoint slope is independent of L, with 
-7 = 9x{L/2) = G"(0.5) « -4.6 (cf. [H); see Fig.Hb), 
where we have also numerically verified ([S]) . We also find 
that the relative sizes of the front and chaotic regions 
remain fixed, with « 5 approximately constant 

and L-independent in the chaotic region (Fig. [Sjc)). For 
these "small" domain sizes, the front translates over long 
times (recall Fig. [5]); interestingly, the statistics of the 
front motion appear consistent with a random walk [ll| , 
as suggested by the trajectories of s{T) shown in Fig. [S] 

The (approximate) scaling form for the time-averaged 
profile g{X) observed for "small" domains breaks down 
for larger L. Instead, for domain sizes in an "intermedi- 
ate" regime with lengths Li < L < L2 ~ 560, the ampli- 
tude of g begins to level off, the front becomes wider and 
less steep, and chaotic fluctuations of A and / decrease 
in amplitude (see Figs. [3H11). Furthermore, the variance 
of the front displacement s{T) decreases strongly with L, 
until the translation becomes imperceptible (Fig. [S]). 

This behavior is transitional to that of "large" do- 
mains L > L2 ~ 560. In this regime the front is sta- 
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FIG. 4: (a) Slope a = gx{0) of averaged profile at midpoint 
of chaotic region. (b) Absolute value of mid-front slope 
7 — —gx{L/2) (□), shown with the mean-square average 
pattern amplitude L'^'^ J^^ {\A{X, dX (0), verifying 
Lengths L and colors are as in Fig. O the vertical lines at 
Li « 220 and L2 ~ 560 indicate approximate transitions 
between "small", "intermediate" and "large" regimes. 
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FIG. 5: Representative trajectories of the front displacement 
s{T) about X = L/2 for the "small" domains L = 38.47r (o) 
and 51.27r (0), and "intermediate" sizes L = 76.87r (a), 
102.47r (□), and 153.67r ([>). 



tionary, s(T) = 0; the amplitude of g{X) saturates at 
max 5 w 62, as does the maximum slope in the chaotic 
region, a = gx{0) ~ 0.5. Indeed, Fig. |3ljb) shows that 
the mean profile g{X) near X = becomes invariant with 
increasing L; this saturation of the profile indicates to us 
that we have reached the large-L asymptotic regime of 
(IT])-(I5]). Since the width of the amplitude death region 
continues to grow with L, while the height is bounded, 
the front slope —7 decays with L, and hence so does 
the amplitude of the fluctuations in A: for large L the 
spatially localized chaotic dynamics superimposed on the 
mean profile are strongly suppressed. 

Transient behavior: The strong L-dependence of the 
properties of the MC equations, within identifiable 
domain size regimes, is apparent also in the tran- 
sient approach to the long-time statistically stationary 
state, as summarized in the time evolution of w{T) — 

[L^^ f{X,T)^ dX]^^'^ — analogous to an interface 
width in the context of surface growth — as in Fig. |6l 
The snapshots from a typical time evolution for a "large" 
domain in Fig. [7] demonstrate an extended coarsening pe- 
riod followed by a remarkable collapse to a single front: 
From small random data, initial growth rapidly estab- 
lishes a sawtooth pattern in /: a concatenation of struc- 
tures, of varying widths and corresponding heights, lo- 
cally reminiscent of the statistically stationary states in 
"small" domains (see Fig. El^a)). Once this metastable 
state of multiple Burgers-like viscous shocks with super- 
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FIG. 6: Evolution of w{T) = [L^^ f{X, Tf dXY''^ for the 
"small" domains L = 25.67r (o) and 51.27r (0), "intermediate" 
domains L — 89.67r (a) and 1287r (<); and "large" domains 
L = 2567r (□) and 307.27r (d>), computed to T = 8 x lO". 
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FIG. 7: (a) Large scale mode / at T = 28000 for the "large" 
domain L = 307. 27r with, for comparison, six copies of the 
L — 51.27r profile from Fig. [T] (b)-(j) Snapshots of / and \ A\ 
(for clarity we plot 5|A|) showing coarsening and collapse to 
a single front for L = 307.27r ^ 964.8. 



imposed chaotic fluctuations is established, a slow coars- 
ening process ensues: front structures grow and merge 
with adjacent fronts, leading to a gradual increase of 
length scales and of w{T) (Fig. [7i;b)-(d)). 

For "small" domains, for which the long-time state has 
the (approximate) scaling form G on average, this coars- 
ening concludes once there is a single front. However, for 
"intermediate" and "large" domains, the gradual growth 
of w{T) through coarsening is followed by a "jump" in 
w{T) (see Fig. IB]) [11], reflecting qualitative changes in 
the profile f{X,T), as seen in Fig. [7Jd)-(h). Specifically, 
having (presumably) exceeded a critical size, one of the 
front structures begins to dominate, and then grows rela- 
tively rapidly by engulfing its neighbors until a state with 
an L-dependent single front is attained. 

Finally, for "large" domains, w{T) overshoots its 
asymptotic value (Fig. [S]), since following the collapse 
to a single front, / is initially nonlinear in the amplitude 
death region, before undergoing slow diffusive relaxation 
(by ([2]) with \A\ =0) to the time-asymptotic linear front 
profile (Fig.[7i;h)-(j)). 

Discussion: The structure of the Matthews-Cox 
equations is reminiscent of that of other well- 

known systems. For instance, viewing ([2]) as a heat equa- 
tion for / with (localized chaotic) forcing, using the heat 
kernel to express / as a quadratic functional of A and 
substituting, the —if A coupling term in Q acts as a 
nonlocal cubic stabilizing term in a GL-type equation. 
Alternatively, in the light of the viscous shock-like behav- 
ior in /, it may be fruitful to view ([2]) as a generalized 
viscous Burgers equation, with a nonlocal forcing term 
determined by ([T]). Such considerations may facilitate 
a theoretical understanding of the unusual behavior we 
have described in the MC equations. 
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